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We consider spatiotemporal chaotic systems^™ for which spatial correlation func- 
tions decay substantially over a length scale £ (the spatial correlation length) that 
is small compared to the system size L. Numerical simulationsUl suggest that such 
systems generally will be extensive, with the fractal dimension D growing in proportion 
to the system volume for sufficiently large systems (L £). Intuitively, extensive 
chaos arises because of spatial disorder. Subsystems that are sufficiently separated in 
space should be uncorrelated and so contribute to the fractal dimension in propor- 
tion to their number. We report here the first numerical calculation that examines 
quantitatively how one important characterization of extensive chaos — the Lyapunov 
dimension density — depends on spatial disorder, as measured by the spatial correla- 

1 



tion length £. Surprisingly, we find that a representative extensively chaotic system 
does not act dynamically as many weakly interacting regions of size £. 

More specifically, researchers have conjectured! that the fractal dimension D of a sufficiently- 
large homogeneous spatiotemporal chaotic system should obey the following simple scaling relation: 

D = C{L/i) d for L » £. (1) 

Here C is a constant and d is the spatial dimensionality of the system, e.g., d = 2 for a large- 
aspect-ratio convection experiment.! Eq. ([l]) follows from an assumption of spatial disorder, that 
parts of a large homogeneous system are uncorrelated and hence dynamically independent when 
separated by distances larger than the length £. The volume L d of the system then acts as a gas of 
N = L d /£ d weakly-interacting regions of volume £ d which contribute independently to the fractal 
dimension. 

Eq. (|H) predicts that spatiotemporal chaos is extensive since the fractal dimension D is pro- 
portional to the system volume L d . This has indeed been verified in numerical calculations of 
one-dimensional chaotic systemsiH and is believed to hold generally. A more subtle prediction 
of Eq. (|]) is that extensive chaos depends on the spatial disorder through a single lengthscale £ 
so that D oc l/£ d for a fixed but sufficiently large system size L. It is useful to formulate this 
second prediction in terms of an intensive parameter so that the system size L and details of the 
boundary conditions no longer enter. This can be done by introducing the intensive dimension den- 
sity S = D/L d , which is the fractal dimension per unit volume for an extensively chaotic system. 
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Eq. (||) then predicts that 

5 oc (2) 

which we test in detail below. Researchers have shown for one case that this quantity is independent 
of the choice of boundary conditions for sufficiently large system sizes (see Fig. 10 of Ref. ||). 

Eqs. (|) and (|) are important to verify and to understand for two reasons. First, they sug- 
gest a simple theoretical explanation for why chaos should be extensive, relating the complexity D 
of a dynamical system (roughly the minimum number of degrees of freedom needed to describe a 
systemtl) to its spatial disorder. Second, Eq. ([[]) suggests a straightforward way to estimate the dy- 
namical complexity of many interesting experimental systems by simply measuring the correlation 
length £.00 This approach would have many benefits compared to time series methods0 since 
large spatiotemporal chaotic systems are high-dimensional, the length of a time series needed to 
estimate fractal dimensions above five is prohibitively large,0 and algorithms for estimating fractal 
dimensions from time series are notoriously unreliable. As an example, we note that scientists 
recently estimated a correlation length of £ ~ 6 mm from voltage measurements on the surface of a 
fibrillating pig heart.0 Since the radius of the heart was R ~ 25 mm, Eq. ([!]) would suggest that 
the fibrillating heart involved approximately D > AttR 2 /tt^ 2 ~ 70 independent regions, a result 
that is consistent with the failure of time-series algorithms to detect low-dimensional chaos. @ 

To investigate the dependence of the dimension density 5 on the correlation length £ we have 



studied chaotic numerical solutions of the one-dimensional complex Ginzburg-Landau equation!!: 

dtu(x, t) = u + (1 + ici)d^u — (1 — ics)\u\ 2 u, (3) 

in a periodic system of size L. The field u(x, t) is complex-valued while the parameters c\ and C3 
are real-valued. Eq. (|3|) is an experimentally relevant0lll continuum dynamical model that holds 
universally near the onset of a Hopf bifurcation from a static homogeneous state to an oscillatory 
state.! It is particularly valuable for testing Eq. (Q) since there is an apparent phase transition in 
large systemsEl — the LI transition — as the parameter C3 is varied for a fixed parameter c\ > 1.8. 
The correlation length £ increases rapidly as one approaches the LI transition from one side and 
so yields a valuable strong variation of £ with parameters. 

We integrated Eq. (|3|) using a pseudospectral method with time-splitting of the operators,!^ 
carefully checking our results for convergence with respect to spatial and temporal resolutions and 
with respect to total integration times. Correlation lengths £ were obtained by examining the 
asymptotic exponential decay Aexp(—x/^) of the magnitude of the complex two-point equal-time 
correlation function C(x) = (u*(x + x' , t)u(x' , t)), where the angle brackets indicate averaging over 
time t and space x'. We calculated a particular fractal dimension, the Lyapunov dimension Dc, via 
the Kaplan- Yorke formula which expresses Dc in terms of the spectrum of Lyapunov exponents.0 
The exponents were obtained by an expensive but widely-used algorithmS that involves integrating 
Eq. (|3|) together with many copies of the linearization of this equation around the solution u. We 
used about 1200 CPU hours on CRAY YMP and on Thinking Machines CM5 supercomputers 
to obtain the results reported below. For most dynamical systems, the Lyapunov dimension is 



believed to be the same as the information dimension D\, which is just one of an infinity of fractal 
dimensions D q that characterize multifractal strange attr actors .0 

Our results are presented in Figures 1-3 for the fixed parameter value c\ = 3.5 and for various 
values of the parameter C3 near the LI transition. We have calculated correlation lengths £ for 
various system sizes up to L = 262, 144 and Lyapunov dimensions Dc for system sizes only up 
to L = 1024 (since the latter is far more expensive to calculate numerically). Figure |](a) confirms 
the prediction of Eq. (|]) concerning extensive chaos by showing that Dc grows linearly with system 
size L over the entire range 64 < L < 1024; for L = 1024, we estimate a large fractal dimension 
of D ~ 152. The dimension density 5 was estimated by a least-squares estimate of the slope of the 
linear portion of this curve. Similar extensive chaos curves were obtained for all C3 values reported 
in this paper. Figure |l|(b) shows how we calculated the correlation length £ from the asymptotic 
exponential decay of the spatial correlation function C(x). Near the LI transition, an unexpectedly 
large system size of L/£ > 1000 was needed to find good exponential decay and to find a value of £ 
that was independent of the system size. 

Fig. [2] summarizes how the Lyapunov dimension density and spatial correlation length vary with 
the parameter value C3 for the particular value c\ = 3.5. As C3 decreases through the approximate 
position of the LI transition (03 ~ 0.74), the Lyapunov dimension density (Fig. ||(a)) decreases 
continuously within our numerical accuracy. The correlation length £ has contrary behavior and 
increases by one order of magnitude close to the LI transition (from 6 to 172). Careful numerical 
tests suggest that £ can not be estimated reliably below C3 < 0.79, probably because it has become 
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too large. As we describe elsewhere in more detail, we find that the £ vs C3 data are well described 
by an algebraic divergence of the form £ oc (03 — c^) a with exponent a = —2.4 and with c\ = 0.735. 

Our key result is given in Fig. [3|, which presents the data of Fig. |2] in a way appropriate for 
testing Eq. @. For spatial dimension d = 1, this equation predicts that the dimension density 5 
should be linear in the quantity l/£ while, in fact, the observed dependence is unambiguously 
nonlinear (the inset figure suggests that a 1/Vf dependence may be more appropriate). If £ is 
truly diverging across the LI transition while 5 changes continuously, Eq. (||) must break down. 
This result is quite surprising since the chaos is extensive (Fig. |l](a)) and an intuitive reason for the 
extensive chaos is that remote parts of the system are uncorrelated and so contribute independently 
to the dynamics. 

We conclude that Eqs. (jlj) and (||) are not correct. One explanation may be that the prefactor C 
in Eq. (|l|) is not constant but depends on £ itself, i.e., details of the local dynamics over a region 
of size £ must be included in the analysis. Another reason may be that there is more than one 
lengthscale needed to characterize Ginzburg-Landau spatiotemporal chaos. A third possibility is 
that Eqs. ^ and (|) are, in fact, correct but involve a new, as-of-yet unidentified, lengthscale A in 
place of the correlation length Further analysis will be needed to distinguish between these 

possibilities and to understand how the dynamical complexity of an extensively chaotic system is 
related to its spatial disorder. 
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Figures 

FIG. 1. (a) Plot of the Lyapunov fractal dimension Dc versus system size L for Eq. (^) with 
periodic boundary conditions. We used the parameters values c\ = 3.5 and C3 = 0.8, a constant 
time step of At = 0.05, a spatial resolution of two Fourier modes per unit length, and a total 
integration time of T = 50, 000 time units. The error in each fractal dimension is smaller than 
the size of the plotted points. 

(b) Log-linear plot of the magnitude of the spatial correlation function C{x) = (u*(x + 
x', t)u(x' , t)) for spatiotemporal chaotic solutions of Eq. @ for the parameter values c\ = 3.5 
and C3 = 0.8. From left to right, the curves correspond to system sizes respectively of L = 1024, 
16384, 32768, 131072, and 262144. A clear exponential decay of correlations is found but only 
for system sizes that are surprisingly large compared to the correlation length of £ = 126. The 
space-time resolution and total integration time were the same as for Part (a). 

FIG. 2. (a) Plot of the Lyapunov dimension density 5 versus the parameter C3 for fixed c\ = 3.5. 
Dimension densities were obtained from the slopes of figures similar to Fig. 0(a). To within 
numerical accuracy, the curve is continuous across the LI transition region. El 

(b) Plot of the spatial correlation length £ versus the parameter C3 for fixed c\ = 3.5. There is an 
apparent divergence which is the signature of long-range order arising from a phase transition. 
The £ values were obtained by least-squares fits to the linear portion of correlation plots similar 
to Fig. 0(b). Systems of increasing sizes up to L = 262, 144 were studied until we had found 
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values of £ that no longer changed upon further increases in L. We were not able to calculate 
correlation lengths reliably for C3 < 0.79 since the value of £ continued to change with increasing 
system size and with increasing integration time. 

FIG. 3. Plot of the Lyapunov dimension density 5 versus inverse correlation length l/£ for the 
parameter value c\ = 3.5. The inset shows the same data plotted versus the inverse square 
root, 1/vT- The dotted lines connecting adjacent points were added only to guide the eye in 
identifying possible linear trends. The nonlinear dependence on l/£ is inconsistent with the 
scaling law Eq. (Q) for a gas of weakly coupled domains that are highly correlated over the 
range of the correlation length £. 
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